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The Faddeev-Yakubowsky equations in configuration space have been solved for the four nucleon 
system. The results with an S-wave interaction model in the isospin approximation are presented. 
They concern the bound and scattering states below the first three-body threshold. The elastic 
phase-shifts for the N+NNN reaction in different (S, T) channels are given and the corresponding 
low energy expansions are discussed. Particular attention is payed to the n+t elastic cross section. 
Its resonant structure is well described in terms of a simple NN interaction. First results concerning 
the S-matrix for the coupled N+NNN-NN+NN channels and the strong deuteron-deuteron scattering 
length are obtained. 
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I. INTRODUCTION 



The four nucleon bound state calculations have, in the last years, reached a high level of accuracy and consistency 
OS , at least as far as the solutions of the corresponding equations are concerned [H-Q . This situation contrasts with the 4N 
scattering problem where despite some pioneering and relevant results |p[-|l4| , there is a manifest lack of convergence 
V/ among the different groups and methods even when using simple interactions. This problem is not only a general 
extension of the 3-body one, in the sustained task of the nuclear Few-Body community to deal with increasingly 
[ complex systems, but we believe it constitutes a qualitative jump in our understanding of nuclear systems. Indeed 
the continuum spectrum of the 4N system (see Figure |l|), with its rich variety of thresholds and structures, provides 
a bridge between the relative simplicity of the A=2,3 problems and the complexity of many-body systems. Even 
when restricted to the energies below the first 3-body breakup threshold, the presence of several resonances at each Z 
channel, the existence of the almost degenerate p+t and n+ 3 He thresholds with, in the middle, the first + excitation 
of the 4 Hc ground state make the understanding of the A=4 chart in terms of fundamental NN interactions an exciting 
and redoubtable theoretical challenge. 

We present here the first solution of the Faddeev-Yakubovsky (FY) equations in configuration space for the four 
nucleon scattering problem. Although the results concerning the bound states ( 4 He and 4 He*) will be discussed in 
some detail, our main interest lies in the 4N continuum spectrum, i.e. the N+NNN elastic scattering and its coupling 
£3 ! to the first inelastic NN+NN threshold. 

The resolution method is based on the angular momentum expansion of the FY amplitudes and the spline expansion 
• i-h , of their radial parts. Orthogonal collocation is used to generate a linear system which is solved by iterative procedures. 
The scattering observables are extracted from a direct inspection of the FY amplitudes in the asymptotic region, in 
H | a natural extension of the methods developed for the three-body case in jl5| . 

The results presented in this paper have all been obtained by using an S-wave NN interaction model and the 
isospin symmetry hypothesis. The Coulomb and mass difference effects are thus not included. This choice, guided by 
methodological reasons, allows a presentation of the formalism and methods in a relatively simple framework. Is is 
remarkable, however, that such a simple model provides a very good description of low energy scattering observables 
even if, as in the n+t case, they are not totally trivial. Some first results including realistic interactions have already 
been reported elsewhere [p| and will be the subject of subsequent publications. 

The paper is organized as follows. In the next section we describe the general formalism and the simplifications 
arising in the case of 4 identical particles. This section contains also the spin, isospin and angular momentum algebra. 
In Section 3 we give some details of the numerical methods used. In Section 4 the results will be presented. They 
include the 4 He ground and first excited state, the elastic phase-shifts and low energy parameters for the N+NNN 
reaction and the N+NNN — > NN+NN first inelastic channel. The energies are restricted below the 3- and consequently 
4-body break-up. Conclusions and perspectives will be given in the last section. 



II. THE FORMALISM 
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A. Faddeev-Yakubovsky equations 



With the aim of solving the Schrodinger equation for N particles interacting via a pairwise potential Vij 

= (1) 



i<j 

Yakubovsky [^7| , generalizing Faddeev's work for N=3 |l^,|l9|], wrote a set of equations whose solutions verify (|l|) 
and which provides a proper mathematical scheme to account for the variety of physical situations involved (see, e.g. 
Figure |^). In the N=4 case, the FY equations can be obtained by first splitting the total wavefunction \& in the usual 
Faddeev amplitudes, ^ij, associated with each interacting pair: 

* = = *12 + *13 + *14 + *23 + *24 + *34, 

and requiring them to be a solution of the system of coupled equations 

{E-Ho)^ =VijJ2^kl (2) 



fc<2 



or equivalently in its integral form 
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(E- 


-H 


- vMj, k 


= Vij 


{E 


-H 




= Vij 


(E- 


-H - 




= Vij 



*y = GoVijV (3) 



G = (E- Ho)- 1 

Each amplitude ^ ij is in its turn split in 3 parts, the FY amplitudes, corresponding to the different asymptotics of 
the remaining two particles: 

*« = *ij,k + + *«,H i < j,k < I (4) 
and obeying the following system of coupled equations : 

+ + (5) 

in which an amplitude ^*>^ 7 , not defined by (Q), has to be understood as being identical to „. Any solution 
of this system of 18 coupled equations, called the FY equations, is a solution of (0) and consequently of the initial 
problem ([!]). Its advantage lies in the possibility to define for system (^|) appropriate boundary conditions ensuring 
the unicity of the solution. Indeed when one of the particle, e.g. labeled by I, is out of reach of the interaction, all 
the amplitudes in (||) tend to zero except k and circular permutations on ijk which obey 

(E-Ho- v^ )fe = Va (¥ ikJ + ¥ jkii ) (6) 

This system of equations, resulting from (||), is equivalent to the 3N Faddeev equations for the particles (ijk). In a 
similar way, when the (ij) and (hi) clusters are free from interaction the only non vanishing amplitudes are ^%j } kt and 
^ki,ij and their corresponding equations in (||) tends to 

{E-Ho- V.j^'.j.u = V«*H,tt (7) 
It is worth noticing that the FY amplitudes can be written in terms of the Faddeev amplitudes in the form: 



= GijVij * w (8) 



where 



2 



@ij — ~ Hq — Vij) 1 
and, according to (|3|), in terms of the total wavefunction: 

*k* = G U Vii Go (V ik + V jh ) * 

= G« VSj Go (9) 

Disregarding the internal degrees of freedom (like spin, isospin), the natural basis for the configuration space is 
provided by the positions of the different particles 

I rir 2 r 3 r4 >=| r \ > <g> | r 2 > <g> \ f 3 > ® | r 4 > (10) 

In order to remove the center of mass motion, it is useful to introduce the relative Jacobi coordinates. Two sets of 
Jacobi coordinates can be defined for each of the 4! arrangements (ijkl). One of K type: 



x K (ijkl) = \l^^-{fj - n) 



y K [i]kl) = \ — —\r k 



m, + m 



, , 2fHjk,l (-. rniTi + rrijfj + m k f k . 

z K {ijkl) = \ J — n J 11 

1 m \ rrii + m,, + m k 



and one of H type: 



x H (ijkl) = \l^^-(rj - ?i) 



m{ijkl) = \l^^-(fi - ffe) 



Mijkl) = 2 J^hi ( m ^ + m ^ _ min + mjfA 
V m \ m k +mi rrii + m,j J 

in which m is an arbitrary mass taken as a reference, and /i Q .^ is the reduced mass of clusters a and (3. However 
some of these 48 coordinate sets are redundant. For instance those obtained by exchanging i ^> j in a K-set or 
i <-> j or/and fe <-» I in the H-set are equivalent. This yields 18 (12K+6H) arbitrary and physically non equivalent 
Jacobi sets, as many as FY amplitudes. Any of these coordinate sets, suitably completed with the center of mass 
coordinate R, constitutes an equivalent description of the 4 particles configuration space. That provides 18 coordinate 
sets and the corresponding bases for the configuration space, equivalent to (|Io|), that will be written in the form 
I XKVKZK^jkVjR > or | XHVHZK{ijkl)R >. The degrees of freedom related to the center of mass motion separate in 
nonrelativistic dynamics and will be hereafter omitted. Although each FY amplitude could be in principle expressed in 
terms of any of these bases, only one of them is appropriate for expanding it. We will denote the resulting components 
by: 



^\j,k{x,y,z) = < x K y K z K (ijkl)\<S> l l:jk > 
*ij,fei(^, y,x) ee < x H y H z H (ijkl)\<^ i:j M > 



The bases described above have to be completed to account for other degrees of freedom like spin, isospin, etc. Further 
details about the formalism and the relation between the different bases sets can be found, e.g., in pofl. 



B. Identical particles 



In the case of four identical particles, the 18 FY amplitudes can be obtained by the action of the transposition 
permutation operators Pij on two of them, arbitrarily chosen provided that one is of K type and the other one of 
H type. We have taken K = ^\ 23 and H = ^12, 34- The 4-body problem is solved by determining the two K, H 
amplitudes which satisfy the following equations: 

(E-H ~ V)K = V [(P 23 + P 13 ) (e + F34) K + e(P 23 + P 13 ) H] (13) 
(E — Hq — V)H = V [(P13P24 + P14P23) K + P13P24 H] (14) 
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where e = ±1 depending on whether the particles are bosons or fermions. The asymptotic equations, i.e. the 
equivalent of (HQ) are in this case 

(E — Hq — V)K = e V(P 23 + P 13 )K (15) 

(E — Hq — V)H = VP 13 P 2i H (16) 

The total wave function is then given by: 

* = * 1+3 + * 2 +2 

* 1+3 = [1 + e(Pia + P 23 )] [1 + e(Pu + P 24 + P34)] K (17) 
* 2 +2 = [1 + e(Pi 3 + P 23 + Pu + P 24 ) + P13P24] H (18) 

Each amplitude $ = K, H is considered as a function of its natural set of Jacobi coordinates x$,y$, z$, defined 
respectively by equations and ( |l2| ) with (ijfcZ) = (1234) and m = mj: 



fx- = r* 2 - ri 




xh = r 2 - n 
Vh = r 4 - r 3 
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They are expanded in angular momentum variables for each coordinate according to 

(j) a (x,y,z) 



< xyz\<& >= 



xyz 



Y a (x,y,z) 



(19) 



where Y Q are generalized tripolar harmonics containing spin, isospin and angular momentum variables and the func- 
tions <p a , the reduced radial FY components, are the unknowns. The label a represents the set of intermediate 
quantum numbers defined in a given coupling scheme and includes the specification for the type of amplitudes (K or 
H). We have used the following couplings, represented in Figure ||: 



K amplitudes : {[(tit 2 ) T J 3 ] T3 U], 



H amplitudes : [(*i*2)t 0> (*3*4)t v ] t ® \ Qx(sis 2 )a al )j x (^( s 3S4)<t v ) ^ 



3*y 



(20) 



where Sj and ti are the spin and isospin of the individual particles and (J V ,T) are respectively the total angular 
momentum, parity, and isospin of the 4-body system. Each component <j) a is thus labeled by a set of 12 quantum 
numbers to which the symmetry properties of the wavefunction impose the additional constraints: (—1)°"^+^+^ — g 
for K and (— l) <7 '»+ r ">+ i » — (—\yy+ T y+ l y — e f or H. The total parity tt is given by in both coupling schemes. 

The radial equations for the components <f> a are obtained by projecting each of equations (|i~^-p~4[) in its natural 
configuration space basis | x$, y$, z$ >. Several steps further |2(J we end with a system of coupled integro-differential 
equations which, most generally, can be written in the form: 



^b aa '4> a '{x,y,z) =^2V aa >(x)^2 f a > a » 4> a "(x f a , a „,y f a , a „,z{ la „) 



+1 



+ / du h a , a »(x,y,z,u) 0«» (a£, a » , y£, a » , ) 



du 



+ 1 



dv g a > a » (x,y,z,u, v) (j) a » (x 9 a , a „ , y a , a „ , z 9 a , a , ) 



(21) 
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with 



h 2 

(E H A Q )<5 QQ ' - V aa '(x) 

m 

„ 2 _ l x (l x + l) 82 _ ly(ly + l) „ 2 _ + 



The functions f a > a ", h a i a », g a > a " contain all the spin, isospin and angular momentum couplings. The arguments 
x a'a" ' x a'a" ' X a' a" >•• • are functions of (x, y, z, u, v) in the more general case, and are detailed in the Appendix. The 
system of equations ( pl| ) has been explicitly written in for the case of four identical bosons. 

The FY components for the different (S,T) channels in the S-wave approximation, i.e. with all orbital angular 
momenta in expansion ( JTgj ) equal to 0, are listed in Table |. In this table, the symbols " — > " and " ~ " denote 
respectively the amplitudes corresponding to an asymptotic N+NNN or NN+NN channel. 

Note that, contrary to the 3N problem, the number of FY components appearing in (|2l]) is infinite even when the 
pair interaction is restricted to a finite number of partial waves. This divergence comes only from the existence of the 
l z additional degree of freedom in the K-like amplitudes. 



C. Boundary conditions 



For all the physical problems we have considered, the boundary conditions can be written in the Dirichlet form. 
The use of reduced radial FY components 4> a in ( |l9| ) imposes for any kind of solution the regularity conditions: 

4> a {x, V, 0) = </> a (x, 0, z) = a (O, y, z) = (22) 

For the bound state problem these conditions are completed by forcing the components 4> a to vanish on the hypercube 
[0,xn] x [0,vn] x [0,zn], i-e: 

(f> a (x,y,Q) = <f> a (x,y,ZN) = (j> a (x,Q,z) = cp a (x,y N ,z) = <j> a (Q,y,z) = 4> a (x N ,y,z) = (23) 

For the scattering problems the boundary conditions are implemented by imposing at large enough values of z the 
asymptotic behavior of the solutions. Thus, for the N+NNN elastic case we impose at zm the solution of the 3N 
problem for all the quantum numbers a a corresponding to the open asymptotic channel a 

(/) aa (x,y,z N ) = t aa (x,y) (24) 

where the functions t aa (x,y) are the Faddeev amplitudes of the 3N problem. Indeed, at large values of z and for 
energies below the first inelastic threshold, the solution of equation ( |l5| ) factorizes into a bound state solution of the 
3N Faddeev equations and a plane wave in the z direction with momentum k a , whereas the solution of ( [L6| ) vanishes. 
One then has, e.g. for an S-wave, 

4> aa (x, y, z) - t Ua (x, y) sm{k a z + 5 a ) (25) 

and the scattering observables are directly extracted from the logarithmic derivative of the K amplitude in the 
asymptotic region: 

k a cot(k a z + 6 a ) = - — ~t~~ rd z (j) aa {x,y,z) (26) 

4>a a {x,y,z) 

Provided we are in the asymptotic domain, the phase shifts thus extracted have to be independent of (x, y, z) and of 
the amplitude index a a - This provides a strong numerical test. An additional advantage of this procedure is that it 
avoids any cumbersome multidimensional integrals. 

In the presence of several open channels, like N+NNN and NN+NN e.g., several resolutions are needed. The 
boundary conditions (|24|) are simply generalized in the form 

<f>a a {x,y,z N ) = X a f aa (x,y) (27) 

in which A Q are arbitrary real numbers and the functions / Qo coincide with t aa if a is a N+NNN channel or are 
analogous to the Faddeev amplitudes for (|l6|) if a is a NN+NN channel. 
Equation ( p4| ) is generalized in the following way: 
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ip aa (x,y, z) = ^-r a h (k a z) + y^ r a/ y "° ^ a ff og ; h + (k a z) J f aa (x,y) (28) 

where /i^ are the regularized Hankel functions and n a is a multiplicity number for the channel a (n a = 4 for 
1+3 and n a = 3 for 2+2 channels) due to the symmetry properties of the total wave function (see, e.g. [p3| ). The 
coefficients r (amplitudes of the incoming waves) and the S-matrix elements are the unknowns. They are obtained 
by identifying with the asymptotic form ( p8| ) the values 0a Q , <9 z (/> Qa at zjq for different solutions corresponding to 
different choices of A's whose number equals the number of channels. We remark that the momenta k a appearing in 
equations (p5|,p6p^) are the conjugate variables of the z-Jacobi coordinates. They are related to the center of mass 



ft k 

kinetic energy E a of channel a according to E a = - — -. The physical momenta, conjugate to the physical intercluster 



distances, are q a = y |fc a or q a = V2k a depending whether a corresponds to a 1+3 or 2+2 channel. 

By the definition ( p8| ) one has unitarity — 1) and symmetry (Sy = Sji) relations. Working with real solutions 
these properties are related (unitarity implies symmetry). However none of them is a trivial consequence of the 
method used but a strong test of numerical accuracy. 

In order that the factorizability takes place the asymptotic 3N or 2N+2N states have to be calculated with the 
same numerical scheme as that used to solve the 4-body problem. This means in practice that they are exact solutions 
of ( |l5| ) or ([l6|), once the z-dynamics are removed. By doing so, the factorization property, valid only in cartesian 
coordinates, is an exact numerical property and leads to stable local results. This behavior is illustrated in Figure |] 
in which the FY amplitudes for an N+NNN elastic scattering at q — and q > are represented. 



III. NUMERICAL METHODS 

The numerical methods used are based on the Hermite spline expansion, orthogonal collocatio n |p^] a nd iterative 
procedures for solving the linear system. An important step in their solution is the tensor trick p|j25|-|27[| . 

We look for the solutions 4> a of the integro-differential system (|2l|) in the form of a tensor product of 1-dimensional 
cubic Hermite splines Si'. 

2Af x + l 27V y + l 2N Z + 1 

<f> a (x,y,z)= ^2 ^2 ^2 c ai j k Si(x)Sj(y)Sk(z) (29) 

i=0 j=0 fc=0 

defined on grids of non-equidistant N q + 1 points G q = {qo, q±, . . . , <7at} where q = x,y, z. A grid G q will be defined 
by giving the number of intervals N q , the end point qjy and the constant scaling factor between two consecutive 
intervals, A q . We will use the following notation G= {N q , A q ,q^} 1 often extended to multi-domain grids according 
to G= {N ql ,A ql ,q N1 ;N q 2,A q2 ,q N2 ; ■ ■ ■} 

The boundary conditions are easily implemented using the properties of the spline functions (value and derivative 
equal to or 1 at the grid points). They result in fixing the values of some of the unknown coefficients c a ijk in the 
expansion (^9|). By exchanging the indices of the two last spline functions in each variable the solution takes the form: 

2AT X 2N V 2N Z 2JV X 2N y 

4> a {x,y, z) = ^2 ^2 ^ c a ijk Sj(x)Sj(y)Sk(z) + f ai jSi(x)Sj(y)S2N :s +i{z) (30) 

i— 1 j — 1 k— 1 i— 1 j — 1 

where f a ij are the coefficients of the asymptotic functions f a (x, y) defined in (p7|). In the cases of a closed channel or 
a bound state these coefficients are zero. 

Collocation points associated with the 3-dimensional grid are used to generate a linear system, the c a ijk being the 
unknowns. The integral terms in ( pi] ) are calculated with a gaussian quadrature rule with typically N u , N v = 6—12 
points in each angular variable u and v. In order to limit the number of parameters we have chosen N u — N v = N g . 

The number of unknowns is given by N = 8N x N y N z N c where N c = Nk + Nh is the number of FY components. 
A rough estimation for the extreme cases of a four bosons bound state and of a scattering state in a realistic problem 
leads to values N c = 2 — 100, N x , N y = 20 — 30, N z = 20 — 40 and consequently to a number of unknowns going 
from N ~ 10 4 toJV~ 10 6 . This implies the use of a dense matrix with ~ 10 12 coefficients. Direct methods are not 
appropriate for solving such huge linear systems, and we have used iterative techniques, which avoid any storage or 
inversion of the matrix. The basic feature of any iterative method is to obtain the solution of the linear system only 
by iterative application of the matrix to an initial vector. This implies a complete calculation of the matrix elements 
each time it is necessary, and requires the intensive use of parallel computers. 
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In the case of scattering states, the boundary conditions fl27|), responsible for the second term in (p0|), generate a 
source term in ( [2l| ) leading to a regular linear system of the type D(E)c = [F + G + H] c + b, where the different 
D, F,G, H matrices are reminiscent of the D aa i , f aa > , g aa i ,h aa i operators. For brevity we will write the system in 
the form Ac = b. 

The numerical method we have chosen to solve this system is GMRES (for Generalized Minimum RESidual algo- 
rithm) HI]. GMRES is a prototype of the so called Krilov subspaces projection methods. Its aim is to minimize the 
residue r — b — Ac of an approximate solution c, starting from a trial vector cb and looking for its correction c — c*o 
in the Krilov subspace JC — {r , Af , A 2 r , • A p ~ 1; Fq} such that the residue r is orthogonal to C = AJC. When the 
dimension p increases, the residue of the approximate solution c can be brought to an arbitrary small value, called 
tolerance. In most of the practical CcLS6S, cl tolerance between 10 3 and 10 6 is sufficient. 

GMRES is a powerful tool when the problem is well conditioned, what is almost never true in a realistic case. The 
way out is to "precondition" the system, i.e. to solve the equivalent problem M~ l Ac = M~ 1 b instead of Ac — b. 
The closer to A is the matrix M, the better is the preconditioning. Our choice was to take the matrix M equal 
to D. As pointed out in p4] , p9| ]l|] its tensor structure, optimized by the choice of cartesian coordinates, allows an 
easy inversion. Our preconditioning technique gives us a converged result after ~ 30 matrix applications, for all the 
considered physical cases and independently of the dimension of the matrix. Examples of convergence curves, i.e. the 
evolution of the residue modulus at each step, are shown in Figure |^. 

In the case of bound state, the asymptotic behavior of the wave function and FY amplitudes leads to a singular 
system D(E)c = [F + G + H]c that will be rewritten for convenience in the form A'c = Ec, the energy E being an 
eigenvalue. It is well known that iterative methods are appropriate for extracting a few eigenvalues, but only those 
with largest modulus. With the 3-dimensional box-like boundary condition (|2^), the existence of overwhelming highly 
positive eigenvalues leaves no hope of obtaining directly the binding energy. 

A crafty trick is to solve the eigenproblem D^ 1 (E) ■ [F + G + H] c = Ac where E is now a parameter ]30| , p^ |. The 
value E will be an eigenvalue of the initial problem when the spectrum {A} contains 1. Furthermore, it can be shown 
by variational considerations that the more excited the state (including the non physical box states), the smaller 
the corresponding A. Thus the eigenvalues of physical interest can be obtained with methods like IRA (Implicit 
Restarted Arnoldi algorithm) |3lj]. We notice that the full inversion of D, including the two-body potential, gives 
a better conditioned spectrum {A} than an inversion of its kinetic term alone, and avoids some of the awkward 
negative eigenvalues generated by the repulsive part of the potential. IRA is a generalization of the power method 
and gives the first eigenvalue in 10 — 15 matrix applications. It is also based on Krilov subspaces, and approximates 
the eigenvalues of A' by those of the restriction of A' to the space spanned by k vectors Xq,x*i, Xf.—i, xq being a 
trial arbitrary vector. Nevertheless, this method requires several calculations for different values of the energy and 
becomes numerically unstable when using highly repulsive two-body potentials, like the interatomic 4 He- 4 He one [ |32[ . 

A more robust technique was finally adopted, often referred as shift-invert method. The initial problem is written 
in the form (A' — Eq)^ 1 ^ — and this technique converges to the energy closest to Eq. It gives very good results with 
a well-balanced mixture of IRA and GMRES. IRA is used to quickly obtain the dominant eigenvalue £o and provides 

the energy E = Eq H . The real difficulty lies in the generation of the corresponding Krilov subspace. It is obviously 

SO 

impossible to invert (A' — E ) since the F,G,H matrices are contained in A'. The step Xk+i = (A' — i? ) _1 a?fc is 
then performed by a GMRES resolution of the equivalent linear system (A' — E )xk+i = Xk- This technique presents 
considerable interest especially for excited states which are easily obtained, independently of the previous convergence 
of the less excited ones. 



IV. RESULTS 



The results presented in this section have been obtained with the spin-dependent S-wave interaction MT I-III: 

= ^expHvr) _ K exp(- Ma r) 



h 2 

The potential parameters and the value — = 41.47 MeV.fm 2 are the same than those used e.g. in [p] j3^j3^1 and are 

slightly different from those given in the original version |50| . Despite its bare simplicity, this potential turns out to 
be very efficient in describing the bulk of low energy properties in the few-nucleon systems. We will first examine 
what we call the S-wave approximation, i.e. the fact that aside from the zero angular momentum of the interaction 
pair, all the angular momenta variables in expansion ( |l9| ) are set equal to zero. The convergence with respect the 
l y , l z expansion will be examined in a second step. 
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A. Bound states 



In the 4N-system, the bound states exist only for the S=T=0 channel. In the S-wave approximation the number 
of FY components is limited to N c = 4 (Nk = Nh = 2) (see Table Q). The binding energies and r.m.s. radius for 
the ground ( 4 He) and first excited ( 4 He*) states are given in Table ||. The corresponding grids are Gi with N g = 12 



and G£ with N g = 6, given in Table III. The estimated accuracy in the binding energies is 0.01 MeV but we notice 
that much less expensive calculations can provide a precise result as well; e.g. the grid G2 with N g = 8 gives also a 
binding of 30.30 MeV. In the ground state we remark a good agreement with the best existing calculations [[l 34|. 

The first excitation which, experimentally, corresponds to a J n = + resonance 0.40 MeV above the p+t threshold 
p5| , appears in this model as a loosely bound state. The binding energy with respect to the N+NNN threshold (E=- 
8.535 MeV in this model) is 0.257 MeV. A similar result was found in Jt| in which different versions of the one term 
separable Yamaguchy potential gives a binding energy varying from 0.07 to 0.40 MeV, depending on their different 
D-state contribution. This + first excitation has been widely considered in the literature as being a breathing-mode 
p(i| p8[ . Our conclusion is however different. We have calculated the regularized two-body correlation functions 
defined by 

Ca.(z)= ]T JJ dydz-\V a ,(x,y,z)\ 2 (31) 

a' (a' x =a x ) 

where ^ a i(x, y, z) represents the total wave function component in the a' quantum numbers, and where a x denotes 
the subset of quantum numbers a relative to the x variable (l x , a x ,j x , t x ). The summation in ( |3l| ) is performed onto 
one of the two basis, K or H. Once the total wavefunction is normalized, the correlation functions are normalized 
according to: 

I dx- C Qx (x) = 1 (32) 

The results are displayed in Figure pi The separated contributions from the singlet and triplet state are plotted 
for (a) triton (b) 4 He ground state and (c) 4 He first excited state. The difference between the correlation function for 
the ground and excited states is remarkable, both in the shape and in the separated singlet-triplet contributions. For 
the excited state one can distinguish the superposition of two structures with two different length scales, the short 
distance part being similar to the triton one. This suggests a 1+3 structure for the 4 He excited state, as can be more 
clearly seen in plot (d) where the results of (c) are compared with those of the triton suitably normalized. Contrary 
to what would happen in a breathing mode, the short distance behavior of the nucleons is that of an unperturbed 
triton with the fourth nucleon being simply a remote spectator. 

By modifying the MT I-III potential strength we failed to pull the state out of the bound region, the 3-body threshold 
moving in the same direction. It seems very difficult for a pure strong interaction to generate a first excitation in 
the continuum. The right position of this resonance is however a crucial point in any attempt to describe the low 
energy data (e.g. p+t) |2Q]. The effect of Coulomb interactions could be enough. However the inclusion of an ad-hoc 
repulsive 4-body term V(p) = Voe~ p can also achieve the same result. 

The preceding results are only slightly modified by the inclusion of higher partial waves in the FY expansion. The 



effects of these contributions can be seen in Table IV. Their smallness shows the validity of the S-wave approximation. 
These results have been obtained using the grid G2 with N g — 8 for the ground state and GJ with N g = 6 for the 
first excitation for which the corresponding triton binding energy is B3 = 8.513 MeV. 



B. Elastic N+NNN scattering 

A crucial point in our method is to ensure the proper description of the 3N asymptotic state. This is used to fix the 
grid parameters for the x, y variables. In order to exhibit the stability of our results, we will compare the phase-shifts 
obtained with several tritons corresponding to increasing numerical accuracies. The considered grids, detailed in 
Table 0, are T 4 with N g = 6 and a binding energy B = 8.593 MeV, T 8 with N g = 8 and B = 8.527 MeV, and Ti 
with N g — 8 and B — 8.535 MeV. We recall that the precise value for the 3N binding energy is B — 8.535 MeV. 
The grid parameters for the z-variable depends substantially on the relative kinetic energy. A zero energy calculation 
requires a relatively large value of zjv but very few points inside are sufficient to describe an asymptotic linear 
behavior. On the contrary as far as the energy increases, the value of zn can be decreased but the oscillations 
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in the relative wave functions demand an increasingly big number of points. A typical grid for the case E = is 
G z = {10, 1.20, 19.0; 03, 1.00, 34.0} 

We have shown in Table [v| the phase-shifts for different (S, T) channels as a function of the center of mass energy. 
For all of them, we have arbitrarily chosen, as in ref. ||, the determination 5{E = 0) = 180°. For the S=0 case the 
comparison between grids T4 and T 8 has been made, showing a good stability despite the fact that grid T4 gives 
only a poor description of the asymptotic state. As it has been already emphasized in section 2, the key point in our 
approach to the scattering problem lies in the coherence between the asymptotic state and the numerical solution of 
the 4N problem, rather than in a very precise description of it. The results corresponding to grid Ts are considered 
as converged. 

A zero energy calculation directly provides us with the scattering length. The results, given in Table VII, show a 
high stability with respect to grid variations (T 8 and T 10 ) and our estimated accuracy is given in column 3. These 
values are in agreement with the existing published calculations Pq| . For T = 1 they are close to those obtained in 
P{ although with the original potential parameters. 



The remaining low energy parameters have been extracted from the phase-shifts and are given in Table VIII for 
the different (S, T) channels . They are defined in the effective range expansion: 



5o(<?) = q ■ cot 6(q) = x(q) 



(33) 



where \(q) = 1 in the usual case but has to be modified in the presence of a near threshold singularity. According to 
p9| , we have taken the form (q real) : 



x(e) 



1 



1- 

\1o 



(34) 



The validity of the expansion (33) in the energy region below the first inelastic threshold is displayed in Figure [?|. 
The usual effective range approximation, i.e. x(q) — 1 and the q A contribution neglected in (|3^), works very well in 
the whole energy domain for all but the (S = 0, T = 0) channel, in which the existence of a near threshold 4 He excited 
state requires the explicit inclusion of the pole contribution ( 34 ) . It is worth noticing that the contributions coming 
from the Vq term are very small and have been included only for completeness. On the contrary the pole contribution, 
existing only in the (S = 0, T = 0) channel, is essential. Expansion ( |3^ ) provides a very accurate parametrization of 
the S-wave scattering amplitude in all the energy region below the break-up threshold. 

We would like to emphasize here the coherence between the 4 He* and the scattering results in the (S = 0, T = 0) 
channel. Inserting the effective range expansion ( |33|) in the S-wave scattering amplitude 



/o(?) 



1 



go(q) - iq 



produces an imaginary pole in the upper complex q-plane with value q* — 0.095i. The corresponding energy E* = 

2 ft * ^ 

= 0.25 MeV is in close agreement with the direct calculation given in Table O. 

3 m , 

A step beyond the S-wave approximation has been taken by keeping the S-wave interaction alone but allowing non 

zero values in the angular momentum expansion of the FY amplitudes. Table IX shows the sensibility of the scattering 
lengths values when the maxima of l y and l z vary from to 2. In the particular interaction model we are considering, 
the number of channels describing the 3N asymptotic state remains the same whereas the number of channels N c of 
the 4N problem considerably increases. These results have been obtained with the grid T 8 and N g = 8 completed 
with a suitable z-grid . The values are well converged with l y , l z = 0, 1 except in the (S = 0, T = 0) case, where the 
big value of the scattering length makes this state rather sensitive to small parameter variations. 



The n+t cross section 



Of particular interest is to apply the preceding results to the description of n+t cross section. On the one hand this 
is a pure T — 1 channel, free from the difficulties related to the Coulomb interactions. On the other hand accurate low 
energy scattering data exist jl0| , showing a structure at neutron laboratory energy Ti a b w 4 MeV which is supposed to 
be created by a series of P-waves resonances |}5| . The calculations discussed in the preceding section have thus been 
completed up to the first 3-body break-up threshold by the inclusion of the first negative parity states J T = 0~, 1~, 2~ 
corresponding to n+t relative P-waves. The resulting total cross sections are plotted in Figure Isl The contributions 







from n+t relative S- and P-waves are distinguished. We notice that the MT I-III model conserves separately L and 
S and consequently the J 7 * = 0~,1~,2~ states coming from L = 1,S = 1 are degenerate. The corresponding cross 
sections differ only by statistical factors. The remaining J w = 1~ state comes from a.n L — 1, S — coupling. In view 
of these results, several remarks are in order: 



1. The scattering lengths obtained in the S-wave approximation (Table [VI Sj ] gives a slightly overestimated value 



for the zero energy cross section <r(0) = 1.85 b. The experimental extrapolated zero energy cross section is 



<r(0) = 1.70 ± 0.03 b [ |4iy4Cf| . However the inclusion of higher partial waves in the FY expansion (Table IX) 
significantly reduces the S-wave approximation result to c(0) = 1.77 b in closer agreement with experiment. 

2. The comparison of the separated spin contribution is not possible since the values of the spin-dependent scat- 



tering lengths, summarized in Table XI, are still controversial 



3. Despite the simplicity of this model, the agreement with experimental data in the resonance region, is very 
good. The n+t P-waves resonances are generated by an NN S-wave interaction alone. An effective 1+3 P-wave 
potential is created due to the exchange mechanism between the four nucleon. We remark however that a first 
attempt to describe this cross section with l z = 1 in the K components but keeping zero all the remaining 
angular momenta in ( [l9| ) failed. The incoming neutron seems thus to be more sensible to the virtual excitations 
of the triton than to the NN pair interactions themselves. 

4. We have calculated the contribution coming from the n+t relative D-waves. They are given by the positive 
parity states L = 2, S = 0, 1. The corresponding phase-shifts at 6 MeV cm. kinetic energy are Ss=o — —3.3° 
and Ss=i — —2.4°, which contributes only a few mb to the total cross section. The results displayed Figure || 
can so be considered as fully converged in the MT I-III model. 



C. First inelastic channel 



The last point to be presented in this work concerns the reaction 

p+ 3 H -> p+ 3 H 
-> n+ 3 He 
-► d+d 

In the isospin approximation employed throughout this work, both n+ 3 He and p+t thresholds are degenerate and 
correspond to the N+NNN of the preceding section. The Pauli principle imposes for the deuteron-deuteron channel 
(— ) i+s+T = +1. If we assume the final state d+d to be in a relative S-wave one has L = T = and J*" = + , 2 + . 
The J 71 " = 2 + state requires a relative angular momentum l z = 2 in the initial N+NNN channel and it is expected to 
give small contribution at very low d+d energy. We will consider only the {J" = 0+,T = 0) state. 

The S-matrix is defined according to (|2^). We present in Table [xjtlie S-matrix elements at several energies, chosen 
with respect to the NN+NN threshold. N+NNN is referred to as channel 1, NN+NN as channel 2. The symmetry 
and unitarity properties are there fulfilled at the level of 10~ 3 , what corresponds to the accuracy of our results. 

With the conventions used above, the scattering cross sections in a given partial wave (J 7r ) in presence of several 
channels a, a', . . . are given by: 

a^ a = ^-^4^-|l-5 aa | 2 (35) 

(J) __ 77 J I a |2 



\Sa'a\ (36) 



where the notation J holds for 2J+1. S ai denotes the spin of the colliding cluster Oj. The corresponding values with 
a, a' = N+NNN, NN+NN are given in Figure ^ (filled circles). We remark no accident in the N+NNN cross section 
when the inelastic NN+NN threshold is open. The NN+NN — > N+NNN cross section displays the usual 1/v law 
of the inelastic process and the crossed channel N+NNN — > NN+NN the expected \f~E law. The elastic d+d cross 
section has been calculated neglecting Coulomb interactions which will make this quantity divergent. However, due 
to the absence of nearthreshold singularity, one can expect small corrections to the low energy parameters obtained 
in this way, like in the n+n versus p+p case. 

Of particular interest is the extraction of the imaginary part of the strong d+d scattering length, a quantity which 
controls the fusion rate in the process d+d — * n+ 3 He. This can be done either from the linear behavior of the d+d 
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phase-shifts J22 = — (art — icti)Q2 (with S22 = e 2 ** 22 ), as displayed in Fig. [lO[ or from the inelasticity in the non 
diagonal S-matrix element which behaves like IS12I 2 = — 4a/(?2 (1 + 2a/g 2 )- Both methods agree with high accuracy 
and lead to the following values: 



a R = +4.91 ±0.02 fm 

a/ = -0.0115 ± 0.0001 fm 



(37) 
(38) 



We remark a very small value of aj that should be only slightly modified once the Coulomb interaction is switched 
on. This small value is due to the small overlapping between the K and H configurations which respectively govern 
the N+NNN and NN+NN asymptotic states. 



We have presented the first solution of Faddeev-Yakubovsky equations in configuration space for the scattering 
states in the 4N system. They concern both the N+NNN elastic scattering and its coupling to the first inelastic 
NN+NN channel. 

The results presented here have been obtained with an S-wave model interaction and in the isospin symmetry 
approximation, i.e. neglecting the Coulomb and mass difference effects. 

The N+NNN elastic phase-shifts have been calculated for the different spin and isospin channels. The low energy 
parameters have been extracted and the validity of the effective range approximation in the energy region below the 
3-body break-up threshold has been analyzed. We have in particular found the importance of including the pole 
contribution in the (S = 0, T = 0) case due to the vicinity of the first = + excitation. We have also found that in 
the framework of our model this + state is bound at 0.25 MeV below the 3N threshold. The study of the two-body 
correlation functions showed that the structure of this state is a 1+3 configuration rather than a breathing mode, as 
it is usually accepted. The coherence between the binding energy of this state and the scattering results has been 
emphasized. 

The n+t scattering cross section has been treated with a special interest and the first negative parity states have 
been included, to account for the structure experimentally observed. The elastic cross section is well described by the 
simple MT I-III interaction, especially in the resonance region. The n+t P-waves resonance is thus reproduced by a 
NN pure S-wave interaction. This shows the difficulty to understand this structure in terms of the NN interaction 
alone. It is created by the direct and exchange mechanism between the incoming neutron and the target nucleons. 

The S-matrix coupling the N+NNN and NN+NN channels has been obtained as well as the corresponding cross 
sections. The analysis of these results allows the extraction of the NN+NN scattering length, whose imaginary part 
controls the fusion process d+d — * n+ 3 He. Its value turned out to be very small (a/ ~ 0.011 fm). 

Further applications of this formalism including Coulomb interactions and more realistic potentials are in process. 
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Grenoblois de Calcul Vectoriel, CEA), and completed on the Cray-T3E at the IDRIS (Institut du Developpement et 
des Ressources en Informatique Scientifique, CNRS). We are grateful to the staff members of these two organizations 
for their kind hospitality and useful advice. 



This appendix aims at the complete expression of the functions /, g, h appearing in ( pl| ) in terms of the corresponding 
quantum numbers a and a' . By projection onto the appropriate K and H basis, the r.h.s. of the equations (11) gives 
raise to integral terms as shown in (13), involving some very complicated functions /, h, g, generated by permutations 
operators. 



The permutation operators Py are completely defined by their action upon each ket of a complete basis. This basis 
is chosen to be an ordered quadruplet |qiq29394 >, where the given i th value represents the state of the i th particle, 



V. CONCLUSION 



APPENDIX: 



1. Permutation operators 



11 



including the space, spin and isospin degrees of freedom. This basis is a generalization of the one given in (|10|). It is 
assumed that corresponds to an exchange between the i th and j set of quantum numbers, e.g. 

P3i\qiq2q3qi >= \q1q2q4q3 > 

More complex permutations are then transparent, e.g. 

•P23-P34I9192Q394 >= \q1q4q2q3 > 

The way the permutation operators act upon the coupled bases (K or H) can be deduced from this basic feature. In 
the configuration space, e.g., a single operator results in a generalized rotation of the Jacobi coordinates. Note that 
because of the symmetry of the bases, the contribution of the operators P13, P13P34 and P14P23 are identical with 
those of P 23 , P23P34 and P13P24. 



2. TL functions 



All the permutation operations can be seen as, at most, two successive rotations, each of them involving only two 
coordinates, e.g. first of all y and z, then x and one of the preceding rotated coordinates. So it is convenient to use 
the functions appearing in the 3-body problem 0], defined as follows. 

Suppose we are to calculate the projection of a given expression [Yj/_ (x')Yi> (ii')]lF(x' , y') 1 F being an arbitrary 
function of x' et y', onto a bipolar-harmonic basis [Yj x (x)Yi where the following relation holds: 



x' \ (a b \ ( x 
y' ={c d \y 



(Al) 



We define TL functions such that: 



dxdy ■ [Y lx (x)Y ly (y)]l ■ [Yl'{x')Yi'{v')]L-F{x',y') = - / du ■ Hf a , ,„ » (x,y,u) ■ F(x',y') 



where, in the integral term, x' et y' are obtained from (Al) with the constraint cos(i,y) = u. These functions are 
given for example by: 

Ht XMX ,(x,y,u)= Yl (_)i'+V+L+fe . a H h c^d X2 i'X'ki Q X 

k,li ,h ,\i ,^2 ,h ,^o(h+h=l' ;Ai+A 2 = A') 



(2l'y.(2X')\ IX ( k l l\ ( k A A \ f h Ai l 



(2i 1 y.(2i 2 ypx 1 y.(2x 2 y. \ o o o^o o o j \o o o 

h A 2 A \ / A I L\j I 1 [ 2 l '\ x^y^ 
){lo A k\[% % ljW PiW 

where I = 21 + 1 et Pj, stands for the A; tft -Legendre polynomial. 

Obviously, in the four-body case, one has to deal with cumbersome recoupling, to isolate bipolar harmonics from 
the more complex coupling scheme (p0|). Nevertheless all the work can be built upon this basic functions. 

3. The only non-zero / function 

fa.a' is non-zero when both a and a' are of H type. Let us define then x^ aa ,, y^ aa i, z^ a , by 





/ Q;Ct / is given by: 

where the isospin contribution t a _ a i is t a ,a' — $t x ,t' 5t' ,t v 1 (— ) Trc+r " _T . 
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4. h functions 



There are two cases where the h functions can be non-zero: when the amplitude a' is of K type. 

a. a of K type 

X 

x aa> ) Vaa' > Z aa> are tncn defined by ( y\ a , I = [ \/I _i o I [ V I and the complete expression for is: 












f £ 




4 1 




S3 




H 


1 ^ 








<7 


J* J 






a 





t_yi+S2-<7x+S2+cr' x -cr J± j s 2 Si (7 X 1 



I Jl n J3 I / ~ f Jl J2 J3 I ( i t -r ~] 

where A < j 4 js je } \ JM.h ■ I H 3b k \, 3 = 2j + 1, t a>a , = (-)*i+*>-t.+**+t.-t 3 . ^ ^ . ^4 ^ ti 

It 1q 7n 7v ?o 7o t 3 3 a; J 



.77 38 39 ) ( 37 38 39 

and with the constraint cos(x,y) = u 



b. a of H type 



x h aa , \ I — ' ■ — \ I r 



aa' 




x aa'i Vaa'i z aa' are now defined by | y\ a , | = | 1 J J y ] and the corresponding expression for 
h aa f is. 



&«,a'(w,«) = *■*»,«'■ E ^^^^-(-) J '-« +, '- J+J '* + ^- J+ '- +, '-'--^( •* H f 2 } 

,, ' l 3y -J 3xy ) 



, i' i' J 1 *\ i ^ Sl ^ I ( 1 1 V 



j'z J J'l J I j£ CT X 4 j I ^ J 2 ^ 

aa' & aa' aa' 

t rp r 3 \ and the constraint cos(x, z) = u. 
5. g functions 

There are two cases where the g functions can be non-zero: the amplitude a must be of K type. 

a. a' of K type 

It is necessary to define in this case x 9 , , y 9 aa i, z 9 aa , and an intermediate coordinate yo such that: 





X aa> 



and | yo 

z 
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The corresponding expression for g aa i is: 



. 1 

g a ,a'{x,y,z,u,v) = - -t 



■<y.l' 






*1 






a' 




h Si 




H 




Si 




L S 






l5 


S 





L z ^ ''yz 



■ A ■ 



l' x A l xy 1 ,_ Yi+s 2 -cr x +s 2 +a'-S A j s 2 «1 0"cc 

S2 < cr [ xyz 
S4 S a' 



9 9 9 
aa' y aoi' ^aa' 



■ n il,v x x(^y^ u )- n xi z z ,i' y i' z (yo^,v) 

with the constraints cos(s, y) = u and cos(yo, z) = v, and where: 

f , - ( _\ti+t2-T x +t 2 +TL-T . a[ *2 t' x T 3 1 ( t 2 h T x 

ta < a ~ [ > A \U T T' z } A \H T 3 t' x 



b. a' of H type 

We define again x 9 aa , , y 9 aa , , z 9 aa , and an intermediate coordinate yo in the following way: 



/ 1 \ / -a \ i -a \ / l Vs 




7a 7§ I ( Vo I and I y I = I £ -I 

\o $ \ ' I \ ' ) \ o o 2 i 




aa' \ / 2 2 U 




The corresponding expression for g aa i is now: 



g a ,a'{x,y,z,u,v) = | • • 





■ A{ l z SA j z ).{-fy+s-,r 3+ L + s-j S l' X y 4 

^_^si+S3— cr' x +a x +S3— a 

xyz 



A 

i'.,. 



So u cr T I ,z 9 , 



• Kl,v x x(^y^)-^xLv y v z (yo,z,v) 

where t Q)a , = (_)*i+*3-^+^+*3-t 3 . A j ^ ^ ^3 | _ ^ | *3 *i ^ | and with thc constraints cos(O) = u and 
cos(j/ ,z) = u. 
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TABLE I. Faddeev-Yakubovsky components for (5, T) 4N 
states in the S-wave approximation. The listed quantum num- 
bers are those defined by the coupling schemes (pp[). The 
symbols '— and '~' emphasize respectively the asymptotic 
N+NNN and NN+NN channels. 



4 He 




4 He[ 


U 


1 


4 He 


4 He* 


3 H 


B 
r.m.s. 


30.31 




30.29 




30.30 
1.44 


8.79 
4.95 


8.53 
1.72 



TABLE II. Binding energies (MeV) and r.m.s. radius (fm) 
for the 4N ground ( 4 He) and first excited ( 4 He*) states. Our 
results for the 4 He binding energy agree very well with the 
best existing calculations. The triton parameters are also 
mentioned for completeness. 



grid Gi 


X 


20 


1.30 


10.0 












y 


15 


1.25 


12.0 












z 


15 


1.25 


15.0 












grid G2 


X 


15 


1.30 


10.0 












y 


10 


1.25 


12.0 












z 


10 


1.25 


15.0 












grid Gt 


X 


08 


1.30 


08.0 


05 


1.10 20.0 








y 


07 


1.30 


10.0 


05 


1.10 30.0 








z 


07 


1.20 


10.0 


13 


1.10 80.0 


10 


1.00 


150.0 



TABLE III. Grids used for 4N ground ( 4 He) and first ex- 
cited ( 4 He*) states. 
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Nk + Nh 




ly 
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B 4 


B% 


B* -B 3 


2 + 2 











30.302 


8.769 


0.257 


8 + 2 





0,1 


0,1 


30.319 


8.763 


0.250 


16 + 3 





0,1,2 


0,1,2 


30.324 


8.770 


0.257 



TABLE IV. Non zero angular momentum contributions to 
4 He and 4 He* binding energies. 
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grid T 8 


: B 3 = 


8.527 MeV 






X 
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1.30 


14.0 
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1.20 
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grid Tio 


: B 3 = 


8.535 MeV 






X 


18 


1.20 


14.0 
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1.00 


18.0 


y 


15 


1.10 


19.0 


04 


1.00 


29.0 



TABLE V. The grids T 4 , T 8 , T i0 used for the tritons 4, 8, 



10. 
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E c (MeV) SJ°) 

S = 0,T = 1 T4 Ts_ 

0.05 169.97 169.95 

0.1 165.85 

0.5 148.99 

1.0 137.13 137.07 

2.0 121.79 

3.0 111.27 

4.0 102.93 102.82 

5.0 96.18 96.13 

6.0 90.18 90.24 

S = 0,T = T4 Ts_ 

0.01 165.17 164.05 

0.02 159.17 157.66 

0.05 147.73 

0.07 142.29 

0.1 135.72 132.98 

0.2 120.52 

0.3 109.52 106.95 

0.5 95.81 

0.7 85.57 

1.0 74.08 71.05 

2.0 49.92 47.16 

3.0 34.71 32.07 

4.0 23.60 21.04 

S = 1,T = 1 T 8 

0.03 172.97 

0.06 170.07 

0.12 166.0 

0.3 158.02 

0.9 142.84 

1.8 129.11 

2.7 119.43 

3.6 111.84 

5.4 100.10 

6.3 94.84 

S = 1,T = Ts_ 

0.01 176.46 

0.02 174.99 

0.05 172.09 

0.1 168.83 

0.5 155.23 

1.0 145.35 

2.0 132.06 

3.0 122.53 

4.0 114.96 

5.0 108.68 

6.0 103.24 



TABLE VI. N+NNN phase-shifts (degrees), as a function 
of the center of mass kinetic energy (MeV), in different (S,T) 
channels. In the S = case, the results with different grids 
are shown. 
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fflgT (fm) 

S T triton 8 triton 10 final value 

I 413 413 413 ± OOT 

1 1 3.73 3.73 3.73 ± 0.01 

14.78 14.76 14.76 ± 0.02 

1 3.25 3.25 3.25 ± 0.01 



TABLE VII. N+NNN scattering lengths values in different 
(S,T) channels, in the S-wave approximation. 
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s 


T 


a (fm) 


r (fm) 


v (fm 3 ) 


go (fm" 1 ) 


a \ 







1 


14.75 


6.75 


0.462 








1 


1 


3.25 


1.82 


0.231 














4.13 


2.01 


0.308 


0.505 


4.0 


1 





3.73 


1.87 


~ 




3.6 



TABLE VIII. Low energy N+NNN parameters, in the 
S-wave approximation. 
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3.25 
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3.08 


17 


3.08 


31 


3.08 


45 





1 


4.13 


4 


4.10 


12 


4.10 


20 


4.10 


28 


1 


1 


3.73 
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3.63 


23 


3.63 


41 


3.63 


59 



TABLE IX. Convergence of low energy N+NNN scattering 
lengths, with respect to increasing internal angular momenta. 



E c cm. 


[Wlc V ) 




S-matrix 






AT i ATATAT 


AT AT i AT AT 
JN IN +1N IN 


O 

Oil 


c 






4.1253 


0.05 


0.772+0.634i 


-0.00473+0.0472i 


l~\ r\r\ a 1-7-1 1 r\ l~\ A TT\ ' 

-0.00471+0.0470i 


0. 882-0. 468i 


4.1753 


0.10 


0.782+0.621i 


0.000488+0.0564i 


0.000486+0.0563i 


0. 771-0. 634i 


4.2753 


0.20 


0.802+0.594i 


0. 0110+0. 0661i 


l~x -1 r\ 1 r\ A/1/11 ■ 

0. 0110+0. 066I1 


0. 566-0. 821i 


4.3253 


0.25 


0.811+0. 581i 


0. 0160+0. 0690i 


0.0161+0. 0691i 


0.472-0.879i 


4.3753 


0.30 


0.820+0.567i 


0.0209+0.0711i 


0.0210+0.0713i 


0.382-0.921i 


4.5753 


0.50 


0.854+0. 514i 


0.0391+0.0748i 


0.0393+0.0752i 


0.0659-0.994i 


5.0753 


1.00 


0.919+0.381i 


0.0759+0.0679i 


0.0758+0.0678i 


-0.480-0.871i 


5.5753 


1.50 


0.961+0.253i 


0.103+0.0508i 


0.102+0.0506i 


-0.784-0.610i 


6.0753 


2.00 


0.983+0.133i 


0.120+0.0292i 


0.123+0.0301i 


-0.934-0.335i 


TABLE X. N+NNN^ 


NN+NN S-matrix elements. The N+NNN channel is labeled by 1, the NN+NN one by 2. 



25 



a ai 
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0.29 


3.13 


± 


0.11 


2.10 


± 


0.31 


4.05 
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0.09 


4.453 
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0.10 


3.325 
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0.016 



[«] 

111 



TABLE XI. Latest experimental results concerning n+t 
singlet (ao) and triplet (ai) scattering lengths (in fm). 
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FIG. 1. The A=4 chart with the more relevant thresholds and resonances (J* ,T). The vertical axis represents a mass scale; 
the horizontal one distinguishes the different values of the electric charge Z. 
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FIG. 2. Different asymptotics to be accounted for in a 1+3 collision. 
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FIG. 3. Spin, isospin and angular momentum coupling schemes used for the K and H Faddeev-Yakubovsky amplitudes. 
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(a) elastic scattering at zero-energy : (f> a (x,y,z) asymptotically linear towards z 

x 




(b) elastic scattering with positive kinetic energy : <f> a (x,y,z) asymptotically oscillating towards z 




FIG. 4. Representations of some asymptotically non-vanishing amplitudes for n+t elastic scattering: (a) zero-energy (b) pos- 
itive kinetic energy. On left, isosurfaces ; on right, sections in a;-variable. The asymptotic factorization between an independent 
pattern on the x, y coordinates and the z-variable motion appears clearly. 
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10 20 30 10 20 30 



dimension of the Krilov subspace 



FIG. 5. Residue modulus obtained with GMRES vs the dimension of the Krilov subspace, i.e. the number of matrix 
applications, for different numbers, n, of unknowns. The linear system Ac — b is normalized such that ||6|| = 1. The initial 
guess is chosen to be 0. 



31 




FIG. 6. Two-body correlation functions C a=c (x) for (a) triton, (b) He ground and (c) first excited states. Solid (dashed) line 
denotes the triplet (singlet) contributions. In (d), the results of 4 He first excited state are compared to the triton correlation 
function suitably scaled. 
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FIG. 7. Effective range expansion, q being the center of mass momentum. The usual one (dashed curve) gives an quite 
accurate description of the scattering amplitude, except in the (S = 0, T = 0) case. The full expansion ( p3| ) (solid line) provides 
a perfect fit of the calculated points. 
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FIG. 8. The n+t cross section a calculated with MT I-III potential is compared to experimental data. The s-wave contribution 
(s, solid line) is slightly overestimated due to the overestimated scattering lengths. The p-wave contribution (p, solid line) 
dominates in the resonance region and is responsible for the very nice agreement (s+p thick line) with the experimental total 
cross section (+'s). The L = 1,S = 1 contribution is split by statistical factors into its J™ = 0~, 1~, 2~ components (dot-dashed, 
dotted, dashed curves), whereas the L = 1,5 = one corresponds to a pure J 7 * = 1~ (dotted curve) partial wave. 
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FIG. 9. Elastic and inelastic cross sections (solid curves), in fm , for the coupled channels N+NNN-NN+NN. They are 
interpolated between the calculated values (filled circles). The energies are given in the center of mass of the incident channel. 
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FIG. 10. Low energy behavior of the deuteron-deuteron phase-shift 8 and determination of the scattering length. Its real 
or imaginary parts an or ai are deduced from the proportionality between the real (solid curve) or imaginary (dashed curve) 
parts of the phase shift and the center of mass momentum qi . 
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